Effect of orbital currents on the restricted optical conductivity sum rule 
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We derive the restricted optical-conductivity sum rule for a model with circulating orbital cur- 
rents. It is shown that an unusual coupling of the vector potential to the interaction term of the 
model Hamiltonian results in a non-standard form of the sum rule. As a consequence, the temper- 
ature dependence of the restricted spectral weight could be compatible with existing experimental 
data for high-T c cuprates above the critical temperature T c . We extend our results to the supercon- 
ducting state, and comment on the differences and analogies between these two symmetry-breaking 
— ■ phenomena. 

-t— > 

Recent results extracted from measurements 0, 0, S [H of the in-plane optical conductivity ct(lj) in 
Bi2Sr2CaCu208+5 (BSCCO) place strong constraints on possible theories of high- T c superconductivity. Special atten- 
tion has been devoted to the partial spectral weight extracted from the optical conductivity in the direction a = x,y 

g; W a (u m ,T) = I Rea aa (cj,T)doj, (1) 

O 



I. INTRODUCTION 



which is analyzed as a function of temperature T and the cutoff frequency u> m , which varies between 1000 cm -1 (0.12 
eV) and 20000 cm -1 (2.5 eV). According to this definition, the weight W includes the condensate peak at oj = 
which develops in the superconducting (SC) state below T c . When ui m is of order of the plasma frequency, up cc 10 4 
\Q • cm -1 , only intraband optical transitions contribute to the measured spectral weight l[T)l. and the so-called restricted 
t**"*- \ or partial sum rule may be applied 0, 0, , which relates W to the average value of the diamagnetic term r aa (see 
CN ! Eq. © below), 

in 

■ W a (cj P ,T) = W(T) = ^ r {T aa ) = 
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where 7ik,cr is the momentum occupation number, V is the unit-cell volume, N is the system size, d = 2 is the 
dimension of the system, e is the electron charge, and we set H = c = 1. The second line of Eq. J2J is obtained under 
the assumption that the interaction term of the Hamiltonian does not couple to the vector potential A, and the final 
equality in Eq. @ is valid only for a nearest-neighbors tight-binding dispersion £k = —2t(cos k x a + cosk y a) with the 
lattice constant a. In this case the spectral weight is a direct measure of the mean kinetic energy K of the system, 

.£h ' an< ^ depends on temperature and interaction strength. For a non-interacting system — /(£k), where £k = £k — M> 
/i is the chemical potential, and f(x) is the Fermi function, so that W(T) increases as the temperature decreases. In 
the presence of a SC instability, the BCS theory predicts that the occupation number is modified below T c and the 

5^ ! partial spectral weight is 



W(T) 1 x 
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(3) 



where Ak = (Ao/2)(cos k x a — cos k y a) is the d-wave SC gap and E^ c — ylk~+ A 2 , is the quasiparticle dispersion 
in the SC state. The kinetic energy increases below T c because of particle-hole mixing, and as a consequence W(T) 
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decreases. The restricted sum rule should be contrasted [9I ITol im IT^ with the full f-sum rule 

„ , . , line 2 , . 

Rea(uj)duj = , (4) 



which relates the integral over all optical transitions (u> m — ► 00) to the total carrier density n and bare mass m, and 
is independent of temperature and interactions. The difference between the full and restricted sum rules is made up 
by transitions between the orbitals described by the low-energy effective tight-binding model and orbitals with the 
energies above cup, not included in this model. As noted in Ref. [loj. there is as yet no complete understanding of the 
relevant orbitals or energy range over which the full sum rule is restored. 

The experimental results and their implications 0,0,0 for the restricted sum rule are the following, (i) Above the 
critical temperature T c , the partial spectral weight W(T) does not decrease when T decreases. Assuming that the 
mean-field Eq. J3J is already valid in the pseudogap state for T > T c , so that W(T) would decrease as T decreases 
even above T c , following Ref. one may conclude that the observed increase of W(T) is in contradiction with the 
opening of a pseudogap. (ii) Using the value t = 0.25 eV, which is typical for cuprates, the tight-binding estimate 
Eq. J2J for the relative thermal variation of W(T) between T c and the room temperature T r , with T r /t ~ 0.1, gives 
(W(T C ) — W(T r ))/W(T r ) ~ 2 • 10~ 3 . This variation appears to be at least one order of magnitude smaller than the 
experimentally observed change in W(T) even for a large u> m = 10 cm -1 . There is even faster increase of W(T) when 
smaller values of ui m are used 0. (iii) Below T c the situation is not clear yet. While early measurements in BSCCO 
samples show that there is an even faster increase of W(T) 0,0, contrary to the prediction of the BCS theory, more 
recent results in BSCCO Q show that there is a flattening of W(T) in underdoped samples for uj m — 8000cm -1 , 
while a BCS behavior below T c is seen in the overdoped BSCCO and in YBCO samples [3|,|5j. 

The possibility of a spectral- weight change below the superconducting critical temperature has been analyzed, for 
example, in Refs. [l3L IliL Il5| in terms of the lowering of the in-plane kinetic energy. In Ref. [l^) the reduction of 
the kinetic energy at T c has been attributed to the transition from a phase-incoherent Cooper pair motion in the 
pseudogap regime above T c (see, e.g. review to a phase coherent motion at T c , while in Ref. 01 a model with 
a frequency dependent scattering rate was used. More recently, the optical conductivity sum rule has been analyzed 
for a model with electron coupled to a single Einstein oscillator ■ 

Here, in contrast to these papers, we focus primarily on the temperature dependence of W(T) above T c and 
on the issue of the compatibility between the pseudogap opening and the absence of a lowering of the spectral 
wei ght . Our purpos e is to show that if the pseudogap originates from a state with circulating orbital currents 
[l8l lig . I20II21I I22I l23L |25L l2g | the opening of the pseudogap can be accompanied by the increase of the partial spectral 
weight. For completeness we extend these results also to the SC state. We find that when a small SC gap opens in 
the presence of a large DDW gap W(T) remains almost constant below T c . 



II. MODEL 



We consider the model with bond currents circulating around elementary plaquettes of copper atoms which is 
described by the effective Hamiltonian 

H = ^][6<cL c k<T + iD k cl a c k+Qa }, (5) 

k.ir 

where c kfT , c ka are creation and annihilation operators for a particle with momentum k and spin a, D k = 
(Do/2)(cosk x a — cosk y a) is the gap, known as the DDW gap j2g, arising from the formation of the state with 
circulating currents, and Q = (7r/a, 7r/a) is the wave vector at which the density-wave ordering takes place. In the 
present paper we do not derive the Hamiltonian by means of a Hartree-Fock analysis of a microscopic model, 
as it has been done elsewhere 0, 0, |2(], 0, 122I l23t l24| . These studies showed that Hubbard-like Hamiltonians 
with additional finite-range repulsion and superexchange interaction can have a stable DDW saddle point. Thus we 
shall parametrize phenomenologically the DDW order parameter and we will analyze within the low-energy effec- 
tive model JSJl the effect of this symmetry breaking on the optical sum rule. It is worth noting that this approach 
has been often adopted in the literature to address several issues related to transport properties in the DDW state 

m m m m m m m m m . 

The derivation of the restricted sum rule depends crucially on the manner in which the vector potential A enters 
the effective low-energy Hamiltonian (0. For lattice models A is usually inserted in a coordinate representation 
by means of the Peierls ansatz 0, 0, Ci — > Cie~ le $ A dr i which modifies the fermionic operator at every site 
i. The dependence of the resulting Hamiltonian on each component A a of the gauge field is H(A a ) w H(0) — 
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J2i e A a {i)ja{i) ~ T^a(*) Tact (*) ' ^ere is the a component of the particle current density and r aa (i) is the 

aa component of the diamagnetic contribution. Thus the total current density j a (i) is j a (i) = —5H/5A a (i) = 
ejP — e 2 T aa (i)A a (i). By evaluating (j a (uj)) in linear response |8|, |3J|, one obtains the complex optical conductivity 



ie 2 



V(uj + 10) 



. 2 

ie 



-((r aa ) -Aaa(w,0)), (6) 



where /Cqq, is the electromagnetic response kernel and the current-current correlation function A QQ (w,q) is defined as 



(7) 



with Sl„ = 2nnT, using the standard analytic continuation iSl„ — > u> + iO. Here j„ (r, q) is the Fourier transform 
of the current density expressed in imaginary-time representation. The Kramers-Kronig relations for the response 
function A aa (cu,0), yield from Eq. (JSJ the optical sum rule (0) for the tight-binding model with nearest-neighbors 
hopping. This derivation of the optical-conductivity sum rule requires the knowledge of r aa , which is easily obtained 
for a Hamiltonian expressed in coordinate representation. For a Hamiltonian in the momentum representation, it is 
more straightforward to apply the sum rule in the form [34| 



1 „ ,. , . P 



Rea(io)dto = — Urn -([p(t,q),^(t,-q)]), (8) 



where p(t, q) is the Fourier transform of the particle density which satisfies the continuity equation 

iq-f(t,q)=0. (9) 



<3p(i,q) , ... 



dt 

We note that substitution in Eq. © of p(t, q) = ^ k a Cj c _ q ; 2 CT Cic+ q /2,<T and of the free-electron expression j p (t, q) = 
(1/to) ^ k a kc k q , 2 .c k+q /2. CT , corresponding to £k = k 2 /2m, returns the full f-sum rule (@J (see |3^). 

In most cases it is assumed that the interaction term of the Hamiltonian involves only density-density coupling, so 
that this is trivially gauge invariant and the Peierls ansatz modifies only the first, "kinetic", term of the Hamiltonian 
(0). For models with nearest-neighbors hopping, r QQ is then related directly to the kinetic energy, (t qq ) = —(K)/d 
and one obtains the usual version of the sum rule given by Eq. pj. However, this assumption is invalid when for 
example "occupation modulated" hopping terms are considered 13] . In particular, if one assumes that the low-energy 
physics of the system can be described by the effective Hamiltonian Jsj, then any distinction in the total energy 
between a kinetic and a potential part is somehow ambiguous. Thus, by transforming the Hamiltonian JSJ to the 
coordinate space, one finds that r aa contains an extra term for Do =/= 0, 

(Taa) = e k(c kCT c k ff ) + «£>k(c k <T C k+Q<r)- (10) 

kcr 

This result is consistent with the derivation JSJ of the sum rule, when one uses the particle current operator compatible 
with the conservation law © and with the equations of motion for the operators c and c\ p^. I29L IsoL l33l| 

fit, q) = K C k-q/2<T C k+q/2 CT ~ *«k C Lq/2 CT C k+Q + q/2 CT ] , (H) 

k,<r 

where i> k = de^/dk and = —dD-^/dk.. The first term of the previous expression relates as usual the particle 
current to the band velocity w k . The second term, which only appears for non- vanishing Do, takes into account 
the contribution of the orbital currents to the electrical conductivity, arising when the DDW order is established. 
Substitution of Eq. JHJ m © yields 

35 - - E wim - (12) 
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FIG. 1: W(T)/t in units of e 2 ira 2 /V [Eq. jnj] for an underdoped (S = 0.13) and an optimally doped (S = 0.16) system. Below 
T c , marked by the arrows, the lower line for each doping represents the spectral weight in the DDW+SC state, Eq. J loft . We 
note that for 5 = 0.13 the decrease of W(T) below T c is almost negligible. Also shown for comparison (dash-dotted line) is 
W(T)/t in the normal state without DDW formation at S = 0.13 [Eq. ©]. Inset: spectral weight plotted as function of (T/t) 2 . 



where £k = yi^ + -D 2 ., and £±.k = —(J-i £k represent the two excitation branches associated with the formation of 
DDW order which breaks translation symmetry. The sum is taken over the reduced Brillouin zone (RBZ). Eq. 112|) 
was derived using the fact that d a v£ = 2ta 2 cosk a a (and d XtV v^ = ±(Do/2)a 2 cos k x , y a), and it reduces to Eq. (0 
for D = 0. 

To extend this result to a state with both DDW and SC order present, we add to the Hamiltonian |5j) an additional d- 
wave mean-field pairing term H p = .[AJJ.c_kiCk| + /i.e.], where Ak = (Aq / 2) (cos k x a — cosk y a). As a consequence, 
the spectral weight in the DDW+SC state reads 



W(T) 
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where -E±.k = 



±,k 



A? is the quasiparticle dispersion in the presence of pairing, and the explicit dependence on 



k has been omitted. Numerical calculation of the spectral weight defined by Eqs. 112|) and (|13[) shows that below 
the temperature Tddw, at which the DDW state is formed, the spectral weight W(T) increases as the temperature 
decreases and the DDW gap opens. As stated above, this increase originates from the second term of Hamiltonian (0 
which effectively enhances the low- frequency conductivity, as already observed in Ref. [^3]. When the temperature 
is lowered further and the SC gap opens at a temperature T c < Tddw, the spectral- weight increase is reduced with 
respect to the DDW state only. 

A more quantitative comparison with experimental data requires the dependences of the gaps Dq and Ao on 
temperature T and doping <5. This issue has been investi gate d within various microscopic models by several authors 
(see for example |23j,|24|). Here, following the lines of Ref. |25|, |27], |32j , and consistently with our effective Hamiltonian 
(|5j), we adopt a mean-field dependence for the DDW gap. We assume [25| that D opens below a doping-dependent 
temperature Tddw(S) = 40[1 — (5/S ) 4 } meV, where 5 is the doping with respect to half filling, and S — 0.2 is the 
critical doping for the DDW formation. We adopt for the temperature dependence of Do the mean-field relation 
Dq(T, 5) — cTD£,w(S)g(T/TDDw(S)), where g(x) = (1 — a; 4 /3)Vl — x 4 , and c is a constant which is used as a fitting 
parameter. To describe the SC transition we solve self-consistently the BCS equations for Ao and /x as functions of 
temperature. We use an analogous set of parameters as in Ref. |25j to estimate the thermal variation of W(T) in 
BSCCO. We focus on an underdoped (5 = 0.13) and an optimally doped (S — 0.16) compound 0,0- The results are 
presented in Fig. 1, where the temperature dependence of W(T) for the tight-binding metal is shown for comparison. 
Below Tddw the weight W(T) increases proportionally to Dq(T), so that the overall increase of W(T) with respect 
to W(Tddw) is more pronounced in the underdoped case (<5 = 0.13), where -Do(0) is larger. In Ref. [lj it has been 
observed that W(T) shows a T 2 dependence above the critical temperature, as expected in the tight-binding model 
Eq. (J2J), but with a much larger slope. To make a comparison with the result of Eq. (T^J, in the inset of Fig. 1 we plot 
W(T) as a function of T 2 . One can see that below Tddw the T 2 temperature dependence of W(T) is still recovered 
over a wide range of temperature, and with a slope in good agreement with the experimental observation. Significant 
deviations are observed approaching T c and below, where our mean-field approach does not reproduce the anomalous 
increasing of W(T) observed in BSCCO in early experiments Q]. Observe however that because at these dopings 
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Ao <C Do, W(T) is not explicitly decreasing below T c , as expected in an ordinary metal-SC transition, but keeps 
almost constant resembling more recent experimental data 0. Finally, we find that the relative variations of W(T) 
below T r = O.lt is of order W{0)/W(T r ) ~ 3 • 1CT 2 , as observed experimentally 0, S 0, and is much larger than 
expected in the simple tight-binding model. 

The previous analysis can be extended to the case where an additional next-nearest neighbor hopping term t' is 
added to the bare band dispersion £k in Eq. (JSJ). Even though both Eq. QlOjl and Eq. H12(l are formally modified, the 
qualitative behavior of the reduced spectral weight is the same, with an increasing of W(T) below the temperature 
for DDW formation. However, as suggested also in Ref. [U, it is likely that the analysis of the DDW state should 
be carried out with a value of t' much smaller then su gges ted by ARPES experiments, leading to small quantitative 
corrections to the previous results. Also, as shown in [l^], the effect of including next-neighbors hopping term t' is 
small, so for qualitative considerations one may consider the model only with a hopping t. 

III. DISCUSSION 

A crucial step in the presented derivation of the sum rule is to use the current operator l|llfl that was considered 
before in Refs. [2?1 12^. l3fl |33|. A different current operator was used instead in Ref. , where it was suggested 
that the gauge field should couple via the Peierls ansatz to the quasiparticle fermionic operators that diagonalize 
the Hamiltonian J^J. As far as the restricted sum rule is concerned, this corresponds to the replacement of the bare 
dispersion law in Eq. @ with the sum of the contributions from the new bands £±(k), and it would produce 
an extra term {2/NV)J2 KBZ (v^D + e) 2 [/(£+) - f(£-)]/E 3 that would be added to the spectral weight (H2J. Its 
contribution to W(T) is negative and of order t y/t/ D$ , as one can check numerically and estimate analytically at low 
doping. The resulting W(T) is then found to decrease below Tddw, m contrast to the experimental observation and 
the result obtained with Eq. I|12|) . Analogously, below T c the ansatz of coupling the gauge field to the quasiparticle 
DDW operators does not reproduced the expression for the superfluid density p s (T) proposed in Refs. |27ll2Sj| . which 
is derived trough the current operator 

The previous discussion shows that there is not yet an agreement in the literature about the proper treatment of 
the transport properties in the DDW state. However, it is worth noting that the form of the current operator and of 
the diamagnetic term used to evaluate the electromagnetic response kernel K, aa (q) in the DDW state are intimately 
related. If a Gauge invariant approximation is used, the response kernel satisfies K, aa {oj = 0, q — > 0) =0 above the 
SC critical temperature T c H,Hj|. This means for example that the diamagnetic contribution (r QQ ) to the superfluid 
density p s (T) cancels the contribution A QQ (if2 ra = 0, q — > 0) providing the vanishing of p s (T) for T > T c . Within the 
low-energy model (JSJ this cancellation holds only if the diamagnetic term l)l(J|l is considered along with the current 
operator (flip, derived from the requirement that the continuity equation 10 be satisfied. The same result does not 
hold by using the mean-field correlation functions defined in Ref. |32(. 

A different approach, which was not investigated here, consists of deriving a proper gauge-invariant approximation 
for the response kernel K- n n b y starting from an underlying microscopic model that provides the basis for the Hamil- 
tonian (0 0, 0, |2jJ, EI HE HE 123 an d including the vertex corrections to the mean-field correlation functions. 
In the case of SC symmetry breaking, one knows that vertex corrections are singular for (w,q) — > 0, satisfying the 
dispersion relations of the collective (phase) mode j^HUl- In the DDW case, where the phase mode is locked by the 
commensurability, vertex corrections are always finite and at zero frequency are related by Ward identities to the k 
derivative of the self-energy associated with the DDW state, i.e. to the term which appears in the definition 
of the current. As a consequence, as observed in Ref. (3j|, the d.c. conductivity, <j(0) calculated with the current 
operator coincides with the exact results for a general, many-body formulation with nonzero vertex corrections. 
This observation suggests that there exists an energy scale below which the approach followed here, where the sum 
rules l|12fl for T > T c and for T <T C were derived directly from the effective Hamiltonian © , describes properly 
the system behavior. However, since there is no straightforward extension of the previous arguments for a(uj ^ 0), 
it is difficult to determine the value of uj m at which the restoration of the more general sum rules © and (QJ should 
be observed in the original microscopic model. In particular, if this cut-off energy resulted to be quite lower than the 
bare plasma edge, the comparison with the experimental data presented before should be reconsidered and referred to 
the data collected up to frequencies lower than 10 cm -1 . As a consequence, the restricted sum rule derived here for 
the mean-field Hamiltonian J5J has yet to be understood from a more general point of view, within a direct analysis of 
a microscopic Hubbard-type model with some short-range interaction that may result in the formation of DDW state. 
As mentioned above, this investigation is rather complicated and cannot be done within the framework considered 
here, so we reserve it for a future work. 

To conclude, we have demonstrated within an effective model that circulating currents can act to modify the 
restricted optical sum rule in a such way that this acquires the same temperature dependence as that observed in 
experiments above T c : the opening of the corresponding gap produces an increase in the spectral weight above T c . 
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Below the SC critical temperature the spectral weight keeps almost constant, as observed recently in Ref. |4J, but 
in contrast with other measurements yj. Since the experimental situation about the behavior of the spectral weight 
in the SC state is not settled, more data are certainly required to definitively establish the possible compatibility 
between our findings and the experiments. As far as our theoretical approach is concerned, we discussed that the 
exact range of validity of this result should still be clarified. Nonetheless, the analysis of the reduced low-energy model 
suggests the possibility that the same kind of deviations from the conventional form of the restricted sum rule could 
be expected in more sophisticated microscopic models (see e.g. |22|). 
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